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The expansion of electromagnetic post-solitons emerging from the interaction of a 30 ps, 3 x 10 18 
W cm -2 laser pulse with an underdense deuterium plasma has been observed up to 100 ps after 
the pulse propagation, when large numbers of post-solitons were seen to remain in the plasma. 
The temporal evolution of the post-solitons has been accurately characterized with a high spatial 
and temporal resolution. The observed expansion is compared to analytical models and three 
dimensional particle-in-cell results providing indication of the polarisation dependence of the post- 
soliton dynamics. 
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PACS numbers: 52.35.Sb, 52.70.Nc, 52.65.Rr 

The propagation of laser pulses through underdense 
plasma is subject to a complex interplay of dispersion 
effects, due to finite electron inertia, and nonlinear it ies, 
due to the relativistic mass increase of the electrons and 
particle redistribution by the laser ponderomotive force 
[l|. This leads to many well-known non-linear phenom- 
ena, including relativistic self-focussing |2| and soliton 
formation [3|46[. 

During its propagation through an underdense plasma, 
the laser experiences a significant energy loss. Being this 
energy loss fully adiabatic, it is mostly translated into a 
red shift of the laser light [7] . In the case of initial plasma 
densities close to the critical density, this frequency de- 
crease may eventually lead the laser to locally experience 
an overcritical plasma, thus being trapped in plasma cav- 
ities. These cavities, whose radius is of the order of the 



electron collisionless skin depth (l e 



c/ujp e , where uj pe 



is the Langmuir plasma frequency) are usually referred 
as electromagnetic (e.m.) sub-cycle solitons [5|- These 
structures tend to be accelerated along plasma density 
gradients [4|, [8[ and, therefore, are slowly propagating, if 
not steady, in a homogeneous plasma. For times longer 
than the ion plasma period, the Coulomb repulsion of the 
ions left inside the cavity causes it to radially expand, 
and the soliton nature is lost: such late-time evolution 
of a soliton is thus commonly referred as a post-soliton 
[9[ . This mechanism is similar to the Coulomb explosion 
that occurs inside laser channels following relativistic self- 
focussing [2[. 

Post-soliton expansion has been analytically studied 
[9|, llO| , using the so-called snowplow model [Uj and the 
isolated spherical resonator model |12J . and numerically, 



using Particle-In-Cell (PIC) codes [9]. Possible experi- 
mental indication of the presence of e.m. soliton struc- 
tures in underdense plasmas has been reported in [13| 
whereas, experimental observation of post-soliton struc- 
tures was first reported in [14j |. where soliton remnants 
were observed in the dense region of a plasma resulting 
from the laser-driven explosion of a thin foil. Due to 
the nature of the plasma employed, clouds of bubble-like 
structures were detected thus preventing a precise charac- 
terisation of the temporal evolution of the post-solitons. 
In this Letter, we report the first experimental observa- 
tion of well isolated post-soliton structures. This allows 
their evolution to be resolved and followed over a sig- 
nificant temporal window. The experimental results are 
compared to analytical and three dimensional (3D) PIC 
code models. 

The experiment was carried out at the Rutherford 




Figure 1: a) Top view of the experimental arrangement, b), 
c) Radiographs of two different shots: outlined with dark 
circles are the density bubbles interpreted to be post-solitons. 



Appleton Laboratory employing the VULCAN Nd:glass 



laser operating in the chirped pulse amplification mode 
[15|. A sketch of the set up is given in Fig. [T] 200 J of 
1 jam laser light contained in a 30 ps Full Width Half 
Maximum pulse ( "Int" in Fig. [1]) were focussed to a peak 
intensity of 3 x 10 18 W cm -2 at the edge of a supersonic 
deuterium gas jet with a backing pressure ranging from 
1 to 100 bar. This resulted, once fully ionized, in an 
electron density of 10 18 -10 20 cm -3 , which is 0.001-0.1 
times the non-relativistic critical density n c . The inter- 
action was diagnosed via the proton radiography tech- 
nique 14j, LL6J], which uses, as a particle probe, a laser 



accelerated proton beam, arising from the interaction of 
a secondary laser pulse (r « 1 ps, E « 100 J, I > 10 19 
W cm -2 , "Bl" in Fig. [1]) with a 20/im thick aluminium 
foil. The virtual point-like source [17j allows imaging 
of the interaction area with a geometrical magnification 
M « (I + L)/l « 6, where I « 6 mm and L « 3 cm, see 
Fig. [TJa. The probe beam, after having passed through 
the gas jet, was recorded by a stack of RadioChromic 
Films (RCF) [3. 

The data set comprised about 30 shots in which both 
the deuterium density and the probe time were varied. 
Two typical proton radiographs obtained at a density of 
0.1n c are shown in Fig. [TJ These images were obtained 
with protons of energy « 4 MeV, 100 ps after the begin- 
ning of the interaction. As a rule of thumb, the electric 
fields are directed from the regions of lighter blue color 
(reduced proton flux) towards the regions of darker blue 
color (increased flux). In both images, a channel created 
by the laser pulse is visible in the low density region at 
the edge of the gas jet. In the dense region inside the 
gas jet a strongly modulated deflection pattern is visi- 
ble along the laser propagation axis. This scaly region 
highlights the presence of a cloud of bubbles that appear 
merged or overlying one another in this 2D projection, 
possibly surrounding the laser-driven channel. Such a 
region visually resembles the cloud of solitons that was 
experimentally and numerically observed in [l4(. Ahead 
of and around this region, isolated bubble-like structures 
are visible (black circles in Fig. [1]), most of them located 
at the end of laser filaments, as numerically predicted 
in [4[. These bubbles are associated with strong probe 
proton depletion with sharp edges. We note that these 
bubble-like structures, which we ascribe to post-solitons, 
were never observed at electron densities of 0.01 n c or less. 

Considering the isolated bubbles allows following the 
fundamental properties of the post-solitons temporal evo- 
lution. Thanks to the multi-frame capability of proton 
radiography [16|, it is possible to follow their temporal 
evolution of these bubbles in the range 80-130 ps after 
the beginning of the interaction, i.e. 40-90 ps after the 
laser has left the gas jet. Both bubbles analyzed were 
found to be effectively stationary in the laboratory refer- 
ence frame and to expand preserving a roughly circular 
shape (Fig. [2]). Bulanov and Pegoraro [10] give analytical 
results for the expansion of ID planar, 2D cylindrical and 
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Figure 2: Zoom in on the bubble structure outlined in Fig. 
[T] in different layers of the RCF stack corresponding to 84 ps 
(a)), 100 ps (b)) and 132 ps (c)). d) Bubble diameter as 
a function of time and fits using results from the snowplow 
model for cylindrical and spherical post-solitons. 



3D spherical post-solitons using the snowplow model. In 
3D the diameter of the sphere is given by do(3 3 ^ 2 t/t s ) 1 ^ 3 
for t ^> t s , where t s is given by ^2-K^noVfiij (Eq), do is 
the initial diameter, no is the initial ion density, rrii is the 
ion mass and (Eq) is the time average of the square of 
the initial oscillating electric field inside the post-soliton. 
We fitted the experimental results with this function tak- 
ing the initial diameter do and the time of creation of the 

1/2 

soliton to as free parameters and (E 2 } = 2 x 10 12 V 
m _1 , which is roughly the average value of 40% of the ini- 
tial laser field. The use of this value for (Eq) is justified 
by recent theoretical results which show that almost 40% 
of the initial laser field is trapped within an e.m. soliton 
[4| . Even though this function was able to fit the experi- 
mental data (see Fig. [2}, it implied a creation time at the 
end of, if not after, the laser pulse duration, which is not 
physically sensible. We therefore tried the 2D cylindrical 
result do(5t/t s ) 2 / 5 : this gave a more physically meaning- 
ful fit with a creation time close to the peak of the laser 
pulse in both cases. The initial diameters from the cylin- 
drical fits are also more reasonable than those from the 
spherical fits, being « 1 jam instead of ~ 3 //m, since we 
have c/ujp e ~ 0.53 jam (see Fig. [2] for the fits result). For 
do ~ 1/im, the time scale of the post-soliton expansion is 
t s ~ 68fs, justifying the assumption t — to ^> t s . Other 
bubbles in Fig. [I] and in different shots (not shown) have 
been found to expand in a similar fashion. 

In order to understand why a 3D structure follows pre- 
dictions for 2 rather than 3 dimensions, we carried out a 
3D run with the PIC code OSIRIS [19]. We considered a 
linearly polarized laser pulse with a wavelength of 1 urn, 



Gaussian spatial and temporal profiles with FWHM of 6 
/im and 1 ps, respectively, and a peak intensity of 3 x 10 18 
W cm -2 incident on a fully ionized deuterium plasma 
with a density of 0.1n c . The interaction was followed for 
5 ps. The simulation box was 350 x 50 x 50 /im, divided 
into 2.4 x 10 8 cells each having 2 particles for electrons 
and 2 for deuterium ions, the time step was 0.196 fs. The 
pulse duration and time considered are 30 times less than 
those in the experiment because of computational limi- 
tations. However, this was long enough for post-soliton 
formation to occur and to follow their expansion after the 
passage of the laser pulse. A number of longer 2D runs, 
with a larger number of particles per cell and a range of 
plasma densities, were also carried out for s-polarization 
(laser electric field out of the plane) and p-polarization 
(laser electric field in the plane). The ion density in Fig. 
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Figure 3: a) Ion density at 3ps from 2D PIC modeling for 
s-polarized light. The dashed lines highlight the post-soliton 
from which the electric field distribution in Fig. [5] was taken. 
b) Iso-surfaces of the ion density at 4.11 ps from 3D PIC mod- 
eling. Sub-channels are formed in the s-plane (perpendicular 
to the laser electric field) following leakage of the laser from 
the channel formed. Several post-solitons are located in these 
sub-channels. 



Ob shows prolate spheroid post-solitons, with an aspect 
ratio of 5:3, lying outside the channel formed by the laser 
in the plane perpendicular to the laser electric field. In 
the 2D runs solitons were only ever formed for s-polarized 
light (see, for instance Fig. Oa). This dependence of soli- 
ton creation upon the laser polarization is in line with 
the PIC code results reported in [4]. The solitons were 
formed as the result of laser leakage due to its breaking- 
up in filaments that led to the creation of sub-channels 
departing from the main channel (as visible in both Figs 
02a and Ob). The post-solitons are elongated along the 
direction of the leaking radiation so that, in principle, 
they can lead to a circular or elliptical 2D projection in 
an experimental arrangement as the one shown in Fig. 
HJa. 

The non-spherical shape of the post-solitons at these 
early times gives a first indication as to why the spherical 
scaling may not apply. To determine the scaling of the 
post-soliton expansion with time, we measured the full 
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Figure 4: Post-soliton diameter X2 extracted from the 3D 
PIC results along with a power law fit and a fit using the 
analytical result for a cylindrical post-soliton. 



width at one quarter of the initial density of the post- 
soliton circled in Fig. [3] along X2 (the smaller diameter) 
at various times, which is given in Fig. HI This defini- 
tion of diameter was chosen because the electric field, 
which determines the shapes seen by the proton probing, 
showed a very high noise level (Fig. [3]). However, the pre- 
cise definition used and the direction in which it is applied 
did not change the scaling. We fitted these points with 
a power law (Fig. HJ, which gave t 0,416 , clearly in better 
agreement with the cylindrical scaling of t 2 l h than the 
spherical scaling of t 1 ^ 3 as the experimental data suggest. 
Following this, we fitted the diameter with the cylindri- 
cal scaling of do[5(t — to)/t s } 2 ^ taking to and do as free 
parameters and (Eq) 1 / 2 = 2.4 x 10 12 V m _1 , taken di- 
rectly from the code results for this post-soliton, which 
is close to the value assumed in fitting the experimental 
results. This gave to = 3.5 ps and do = 0.53 /im, which 
gives t s = 31 fs (Fig. [4j). Here the assumption of the 
analytical model that t — to » t s is again clearly satis- 
fied. The laser was still within the computational box at 
3.5 ps. The value of do appears somewhat low, since it 
is less than 2/ e , but this is due to the way the diameter 
was determined, which clearly gave a value smaller than 
the outer edge of the post-soliton, which was difficult to 
determine unambiguously from the 3D results. 

The observed polarization dependence of soliton for- 
mation suggests why the spherical scaling does not ap- 
ply. The snowplow model assumes in fact total reflection 
of the trapped light at the overcritical soliton wall and 
therefore no plasma heating. This is a good approxima- 
tion only for s-polarized light, since a p-polarized wave 
will be indeed absorbed by the overcritical walls [20|. 
This might explain why p-polarized light is not able to 
excite stable e.m. solitons [4|. 

The electric field structure in the post-soliton was 
much clearer in the 2D PIC results. Fig. [5] shows a cross- 



cut through the electric field of the post-soliton outlined 
by the dashed lines in Fig. 03b, and the corresponding 
charge density. The electric field shows two sharp, bipo- 
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Figure 5: a) Cross-cut of the post-soliton charge density from 
2D PIC modelling at 3 ps. b) Corresponding electric field 
distribution, c) Comparison between the modulation of the 
proton probing beam, as directly obtained from the RCF (red 
line) and the modulation simulated via the PT code (blue 
line) . 



experimental results is shown in Fig. [5j the two curves 
show good quantitative agreement, further confirmation 
of the interpretation of the structures as post-solitons. 

In summary, we have reported the first experimental 
measurements of the late time expansion of post-solitons 
following the propagation of a relativistically intense laser 
pulse through an underdense plasma. The post-soliton 
expansion has been temporally resolved with high tem- 
poral and spatial resolution. The post-soliton expansion 
is best described by the analytical prediction for cylin- 
drical, not spherical, post-solitons. 3D PIC code results 
show the same behavior and, consistently with 2D re- 
sults, indicate that this is due to polarization effects. 
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tory project and was supported by the UK Science and 
Technology Facilities Council. 



lar peaks located at the edge of the post-soliton with 
maximum amplitude, at 3ps after the beginning of inter- 
action, of 2 x 10 11 V m _1 ; this field would deflect the 
probe protons in the experiment. In order to compare 
this result with the experimental ones, we used a particle 
tracing (PT) code that traces proton trajectories from 
a point-like source through a given time-dependent 3D 
electric field distribution up to the proton detector, giv- 
ing a 2D proton density map at the detector plane. The 
electric field profile used as an input for these simulations 
was the one produced by the 2D run with a spherically 
symmetric distribution around its center. The electro- 
static field is expected to scale as E 2 ~ d~ 4 ; this is con- 
cluded on the basis that the ratio between the energy and 
the frequency of the trapped electric field is an adiabatic 
invariant during the post-soliton expansion [9]. Given a 
cylindrical scaling for the post-soliton diameter, we ex- 
pect thus the amplitude of the electric field to decrease in 
time as E ~ £ -4 / 5 . At 100 ps, the time of the measure- 
ments, the maximum amplitude of the electrostatic field 
should thus be 2 x 10 n (i00/3)- 4 / 5 « 10 10 V m" 1 . Such 
a strong electric field would induce an almost complete 
depletion of the probing protons in correspondence to the 
post-soliton structure, as the experimental raw data in- 
dicate. Therefore, the initial conditions for the PT code 
were a uniform proton beam with an energy of 4 MeV 
crossing an electrostatic field with a spatial profile as the 
one in Fig. [5j with spherical symmetry and a maximum 
amplitude of 10 10 V m _1 ; the bipolar peaks of the elec- 
tric field extended for 10/im each, separated by a plateau 
region 60/im long. A cross-cut of the code results and the 
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